Late-time attractor for the cubic nonlinear wave equation 
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We apply our recently developed scaling technique for obtaining late-time asymptotics 
to the cubic nonlinear wave equation and explain appearance and approach to the two- 
parameter attractor found recently by Bizon and Zenginoglu. 

I In [l| Bizon and Zenginoglu conjecture and present some analytical and numerical evidence that 

the spherically symmetric cubic nonlinear wave equation in three dimensions 



au = u^ (a^d^-A) (1) 

. has a universal two-parameter attractor for the late-time asymptotics 



Ph. 



for a big family of initial data, being itself an exact solution of ([T|). The rate of approach to the 
attractor ([2]) with suitably chosen a, b is for a fixed r. The aim of this Article is to prove 
this conjecture by deriving aprecise late-time asymptotics and comparing it with the form of the 
attractor ([2]). Recently, in [4], we have developed a method for obtaining late-time asymptotics 
from the scaling properties of a given wave equation. It allows us to show that, at least for small 
I initial data, the generic late-time asymptotics of solutions to ([T]) takes the form 

u ^ u ^ , , Alt , A2i3t^+r^) , 

^: ^(^'^) = ^"o(t,r) + ^2372 + (^2372)2 + -(^^372)^ + ••• 



^ where Ai are given in terms of the initial data. This expansion coincides with the attractor 

Q ■ when Aq = V^/b and Ai = \/2/(2a6 -|- 1) up to the term of order — r^)~^. Hence, as long 

On . as Aq,Ai are non-zero there exist unique parameters a,b such that the solution u asymptotically 

' approaches the two-parameter attractor ([2]). The rate of approach is determined by the failure of 

([2]) to reproduce the third term in the generic asymptotics ([3]) which contains a third independent 
^ ■ parameter A2. 

Such precise asymptotic analysis has only become possible with the recent developments in 
the perturbation theory. Crucial are the first rigorous works [1, 0] allowing for establishing a link 
between the late-time asymptotics and the small initial data (see also references therein for earlier 
non-rigorous but important works, e.g. by Bizon et al). They build up on earlier decay estimates 
for small data by John [sl and Asakura [a]. Later, in [2], an equivalent technique based on scaling 
has been introduced which simplifies the asymptotic calculations of higher order terms. In it, the 
initial value problem for a class of nonlinear wave equations is considered which we restrict here 
for the sake of simplicity to 

□u = uP (4) 

with integer p > 3 and small initial data 

u(0,r) = e/(r), dtu{0,r) = eg{r) (5) 

being smooth functions of compact support in three spatial dimensions restricted to spherical 
symmetry. Then, for small e, we have ueC°°{R^+^). 
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The main idea of deriving asymptotics from scaling is based on the observation of Lindblad 0] 
that in the Umit e — t- the solutions Us of (jl])-© tend, under suitable scaling, to some nontrivial 
It* which satisfies a linear wave equation with a distributional source. This equation can be solved 
exactly. For small but finite values of s the solutions Ue are near to in a suitable sense with a 
uniform error bound such that determines the late-time asymptotics of u^. 

A straightforward generalization of Theorem 1 of [2] leads to 



n 



r 

k=0 



1 1 



+ o{- . +0 



{t + r){t- r)P+"-2 J y(^t + r){t- r)P 



0-2 



(6) 



for e — )• 0, t — r > a given scaling parameter < a < p{p — l)/(p + 1) and any nonnegative 

integer n. uq solves the corresponding linear problem ©-([H]) (see below) and ^ are determined 
by the initial data and are defined below. 

The error terms mean that the asymptotics holds w.r.t. to weighted-L°° norms (cf. [2]). Here, 
they are restricted to the region t — r > l/e"" and imply a uniform convergence there, as e — )• 0. 
The first error term describes correction entering at the same nonlinear order (the same power 
of e) as the leading terms but having faster decay in time while the second error term, with 
^0 '■= {p — 1)(1 — a) + a[{p — 1)^ — 2]/p, stays for corrections with the same decay in time as the 
leading terms but entering with higher powers of e. 

For p = 3 and n = 2 we essentially obtain the asymptotic expansion ([3]). The only point 
which requires some additional but straightforward work is to show that all expressions appearing 
in the second error term are actually of the same functional form as those alredy present in the 
leading asymptotics multiplied by higher powers of e. Then they don't change the character of the 
asymptotic expansion but only alter the constants, thus leading to ([3]) with Ai^{e) = B^^j^e^ + 0{e'^). 

In [ll, Sec. 5.1.2] also non-generic solutions with faster late-time decay than this given by ([2]) 
have been found numerically. This can be explained by the observation that the initial data can 
be chosen such that at late times Aq = or Aq = Ai = 0, etc., thus leading to faster decay in ([3]). 



The scaling technique 

Here, we briefly introduce the method of scaling developed in and extend it to calculate higher 
order terms appearing in ([6]). All technical details and proofs can be found or easily adapted from 
0]. The scaling method has an advantage over the standard perturbation theory in generating 
simpler-to-solve effective equations for the higher order asymptotic terms. 

In the first step we solve a corresponding linear equation with removed scale factor e 

□no = 0, (7) 

no(0,r) = /(r), dtUo{0,r) = g{r). (8) 

Its solution can be written in the form 

, , hit - r) - h{t + r) 
Uo{t, r) = ^ ^-—^ (9) 



where 

h(T.'] ■= - , 

2-" ' 2 



h{x) := -^f{x) - \ I yg{y)dy (10) 
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has compact support (the functions f{r),g{r) have been symmetrically continued to negative r). 
Next, we subtract the linear solution from the nonlinear one by introducing 



which satisfies 



w{t, r) := u{t, r) — euo{t, r) 



Now, we scale this function to 

We{t,r) := e-^wie-^t^e-^r) 
with b = p + a{p — 1) and some a > to be chosen later. It satisfies 

nWe{t,r) = e-" \e^'P-^^^''^v-^'^We{t,r) + e-''uo{e-''t,e-''r) 



(11) 
(12) 

(13) 



(14) 



For this equation we want to consider the limit e — )• 0. 

Let us recall the following fact of the distributional calculus: any smooth function H{x) of 
compact support can be squeezed to the delta distribution under appropriate scaling as e — )• 0. 
The corrections can be written as a sum over derivatives of the delta 



e-^H{e-^x) ^ Cq5{x) + eCi5'{x) + . . . 
what, in the precise (distributional) sense, means 



lim 

e-i>0 



/I 1 

k=0 



where g G is a test function, n is any nonnegative integer and '■= f x^H{x) dx. 

Having this in mind we observe that the first term in (jl4p . by use of the representation 
have a distributional limit (in the above sense) 

£ [e UQ{e t,e r)\ ~ Oq,o — he : 1- • • • 



where 



C, 



q,i 



x'h'^{x)dx. 



(15) 



(16) 



will 



(17) 



(18) 



The terms S'^^'^ {t + r) will further play no role since their support is outside of the region of our 
interest t + r > 0. 

The sum in p4p will be treated as an error term. The expansion 



We = W^^^ + e"VF(^) + • • • + e'^'^VF^") + e("+i)«^ + We 
allows us to write the following limiting equations 

(5('=)(t-r) 



rP 



(19) 



(20) 
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for k = 0, ...,n which can be solved exactly 

e(t - r) 



B, 



p,k' 



-2 



(21) 



with Bp^k ■■= 2P+^'^Cp^k/{p + k- 2). 

The first error term W comes from the truncation of the expansion ()17p of 
[s~"'UQ{e~"'t, s~"'r)y into deltas and can be bound by 



{t + r){t-r)P+''-^\W{t,r)\ < C 



while the second error term satisfies 



k=i ^ ^ 

In [2] we have shown that it can be bound 

+ r){t- ry-^\We{t, r)\=0 (^e^P-m-a)+a[{p-ir^2]/p^^ 

uniformly for all t — r > 1. It vanishes in the limit e — )• for < a < p{p — 1) /{p + 1) 
It allows to write the asymptotics of u 

u{t, r) = euo{t, r) + eP+'''^P-'^'>We{e% eV) 

G(t-r) 



euo{t,r) + 



'p,k- 



k=0 



+ 



{t + r){t-r)P+'^ 



-2 



+ 



-P+Ao 



{t + r){t-r)P~ 



(22) 



(23) 



(24) 



(25) 



which holds uniformly in t — r > e~"' for e — )• 0, a given scaling parameter < a < p{p — l)/(p + 1) 
any nonnegative integer n and Aq := (p — 1)(1 — a) + a[{p — 1)^ — 2]/p. 
This procedure can be made rigorous as has been demonstrated in 
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